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Abstract. 

Numerical simulations for the merger of binary neutron stars are performed in 
full general relativity incorporating both nucleonic and hyperonic finite-temperature 
equations of state (EOS) and neutrino cooling. It is found that for the nucleonic 
and hyperonic EOS, a hyper massive neutron star (HMNS) with a long lifetime 
(^lifo Jc^ 10 ms) is the outcome for the total mass w 2.7 Af©. For the total mass 
w 3 Mq, a long-lived (short-lived with « 3 ms) HMNS is the outcome for the 
nucleonic (hyperonic) EOS. It is shown that the typical total neutrino luminosity of 
the HMNS is ~ 3 - 6 x 10'^'^ erg/s and the effective amplitude of gravitational waves 
from the HMNS is 1 - 4 x IQ-^a at /« 2 - 3.2 kHz for a source of distance of 100 
Mpc. During the HMNS phase, characteristic frequencies of gravitational waves shift 
to a higher frequency for the hyperonic EOS in contrast to the nucleonic EOS in which 
they remain constant approximately. Our finding suggests that the effects of hyperons 
are well imprinted in gravitational wave and its detection will give us a potential 
opportunity to explore the composition of the neutron star matter. We present the 
neutrino luminosity curve when a black hole is formed as well. 
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1. Introduction 

Coalescence of binary neutron stars (BNS) is drawing attention of researchers in various 
fields because it is one of most promising source for next generation kilo-meter-size 
gravitational-wave (GW) detectors [1] and a possible candidate for the progenitor 
of short-hard gamma-ray bursts [2] as well as a high-end laboratory of the nuclear 
theory [3]. 

BNSs evolve due to gravitational radiation reaction and eventually merge. Before 
the merger sets in, each neutron star (NS) is cold, (i.e., thermal energy of constituent 
nucleons is much smaller than the Fermi energy), because thermal energy inside the NSs 
is significantly reduced by neutrino and photon coolings due to the long-term evolution 
(typically > 10^ yrs) until the merger. By contrast, after the merger, shocks are 
generated by hydrodynamic interactions, which heat a merged object, hyper massive 
neutron star (HMNS), up to ~ 30 — 50 MeV, and hence, copious neutrinos are emitted. 
The rest mass density p of the merged object will go well beyond the nuclear matter 
density pnuc ~ 2.8 x lO^^'g/cm^ after the merger. In such a high density region, properties 
of the NS matter are still poorly understood. In recent years, the influence of non- 
nucleonic degrees of freedom, such as hyperons, meson condensations, and quarks, on 
properties of the NS matter has been discussed extensively. Among exotic particles, 
A hyperon are believed to appear first in (cold) NS around the rest mass density of 
p ~ 2 — 3 Pnuc [31 H] . The presence of such the exotic particle results in softening of 
equation of state (EOS). Hence, if it exists or appears in the course of BNS mergers, it 
will affect the dynamics of BNS mergers and be imprinted in emitting GWs as a result. 

Motivated by these facts, numerical simulations have been extensively performed 
for the merger of BNS in the framework of full general relativity in the past decade since 
the first success in 2000 [5]. In particular, it is mandatory to do numerical relativity 
simulations of BNS mergers implementing a finite temperature EOS as well as a neutrino 
cooling. In this proceeding, we present the results based on Refs. [6l [7], in which the 
NSs are modeled with a nuclenoic or hyperonic EOS and the cooling is properly taken 
into account. Specifically, we adopt Shen-EOS for the nucleonic EOS [8] and a finite- 
temperature EOS including contributions of A hyperons as the hyperonic EOS [9j. In 
the following, we report a comparison between the models with these EOSs, particularly 
focused on the dynamics of BNS merger, emitted GWs, and neutrinos. 

2. Framework 

2.1. Formulation and EOS 

Numerical simulations in full general relativity are performed using the following 
formulation and numerical schemes: Einstein's equations are solved in the 
original version of Baumgarte-Shapiro-Shibata-Nakamura formulation implementing the 
puncture method |T0]; a fourth-order finite differencing in space and a fourth-order 
Runge-Kutta time integration are used; a conservative shock capturing scheme with 
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third-order accuracy in space and fourth-order accuracy in time is employed for solving 
hydrodynamic equations |6l[7]. We solve evolution equations for neutrino (Y^), electron 
(Ye), and total lepton (Yi) fractions per baryon, taking into account weak interaction 
processes and neutrino cooling employing a general relativistic leakage scheme for 
electron type (ue), electron anti-type (ug), and other types (fi/r) of neutrinos (ux) 
[TT| IT2] . We implement Shen-EOS [8], tabulated in terms of the rest-mass density 
(p), temperature (T), and Y^. or YJ, as the nucleonic EOS and a finite-temperature EOS 
including contributions of A hyperons (Hyp-EOS) as the hyperonic EOS [9]. These EOSs 
produce a maximum mass of zero-temperature spherical NS of ~ 2.2Mq for Shen-EOS 
and ~ I.SMq for Hyp-EOS. Although the maximum mass of Hyp-EOS is slightly smaller 
than the observational constraint of PSR J1614-2230 (Mji6i4-223o = 1.97±O.O4M0) |13], 
it deserves employing to explore the impact of hyperons on the BNS merger because 
it can be a viable of EOS of the NS matter not in extremely high densities and hence 
for studying the evolution of the HMNS that do not have the extremely high densities 
for most of their lifetime. It should be also noted that it is possible to produce a more 
massive neutron star with hyperons by fine-tuning the parameters (see e.g, in Ref. jl4j). 

2.2. Numerical issues and set up 

Numerical simulations are performed preparing a non-uniform grid as in Ref. P, [TJ [IHl 
[16] . The inner domain is composed of a finer uniform grid and the outer domain of 
a coarser nonuniform grid. The grid resolution in the inner zone is typically chosen 
to cover the major diameter of each NS in the inspiral orbit by 80 grid points. Outer 
boundaries are located in a local wave zone at ~ 560 — 600km, which is longer than 
gravitational wavelengths in the inspiral phase. As a monitor of the accuracy, we 
check the conservation of baryon rest-mass, total gravitational mass, and total angular 
momentum and find that they are preserved within the accuracy of 0.5%, 1%, and 3%, 
respectively. 

References P, [7] focus on the merger of equal mass BNS with total mass ranging 
from 2.7 Mq to 3.2 Mq. In this proceeding, we report the models with total mass 
^2.7 Mq and ~ 3 Mq. There are two possible fates of BNS: If its total mass M is 
greater than a critical mass Mc, a black hole (BH) will be formed soon after the onset 
of the merger, while a differentially rotating HMNS will be formed for M < M^. M^ 
for Shen-EOS (Hyp-EOS) is 2.8 - 2.9 Mq (2.3 - 2.4 Mq). Therefore, for the binaries of 
M 2.7 Mq, we expect that Shen-EOS model will form a HMNS and Hyp-EOS model 
will collapse to a BH. The models with M f» 3 Mq are anticipated to collapse to the 
BH, irrespective of the EOSs. We name Shen-EOS (Hyp-EOS) model with M ^ 2.7 Mq 
S135 (H135) and those with M ^ 3 M© S15 (H15). It should be noted that these figures 
of the mass are motivated by the narrow mass distribution of observed binary neutron 
stars. 
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3. Numerical results 

Figure [Uplots the evolution of maximum rest-mass density, pmax, maximum temperature, 
Tmax, and maximum mass fraction of hyperons, XA,max as functions of t — tmerge where 
Emerge is au approximate onset time of the merger. Before the merger (t < tmerge), Pmax 
and Tjnax for H135 (H15) agree well with those for the corresponding S135 (S15), because 
-^A,max in this phase is small as 0(10^^) and effects of hyperons on dynamics are not 
significant. After the merger sets in (t > tmerge), on the other hand, XA,max increases to 
be > 0.1 in accordance with the increase in pmax, and hyperons play a substantial role 
in the post-merger dynamics. The panel (b) in Fig. [1] tells us that Tmax in the HMNSs 
reach to 50 - 60 MeV just after the first contact, irrespective of the models. For H15, it 
rapidly increases to a high value of 130 - 140 MeV just before the collapse to a BH. For 
H135 and Shen-EOS models, it decreases due to the neutrino cooling in the subsequent 
evolution and then relaxes to ~ 25 - 40 MeV. We observe the rapid enhancement of 
^max just before a BH formation in H135. 

Although the total mass, M, is larger than the maximum mass of the zero- 
temperature spherical NS for all the models, a HMNS is formed after the merger, 
supported by the centrifugal force and thermal contribution to the pressure [HI IT7] . 
In particular, the enhancement of the critical mass due to the thermal effect is 
prominent in S15, which does not collapse to a BH during the simulation contrary to 
our expectation. The HMNSs subsequently contract by emission of GWs, which carry 
energy and angular momentum from the HMNS; pmax increases in the gravitational 
radiation timescale. For Hyp-EOS models, they collapse to the BH at t = ^bh where 
^BH — ^mcrgc ~ H-O ms for H135 and ~ 3.1 ms for H15. After the BH formation of Hyp- 
EOS models, an accretion torus is formed around the BH and gradually settles down 
to a quasi-steady state. The torus mass is ~ 0.082, and 0.035 Mq for H135 and H15, 
respectively [7]. The torus mass of H135 is greater than that of H15 due to the longer 
transient HMNS phase, in which the angular momentum is transported outwardly. Note 
that the lifetime of HMNS for H135 is longer than that for H15. At t — tmerge ~ 20 ms 
for Shen-EOS models, the degree of its nonaxial symmetry becomes low enough that 
the emissivity of GWs is significantly reduced. Because no dissipation process except 
for the neutrino cooling is present, the HMNS has not yet collapsed to a BH. We expect 
it will be alive at least for a cooling time, tcooi = Et^/L^ ~ 2-3 s, where -Eth is total 
thermal energy and is total neutrino luminosity [6]. 

Figure [2] plots neutrino luminosities as functions of time for three flavors {ug, Pg, 
and sum of Ux). It is found that electron anti neutrinos are dominantly emitted for any 
model. The reason for this is as follows: The HMNS has a high temperature, and hence, 
electron-positron pairs are efficiently produced from thermal photons, in particular in 
its envelope. The positron capture n + ^ p + Ue proceeds more preferentially than 
the electron capture p + e~ n + because the proton fractions is much smaller 
than the neutron fraction. Such hierarchy in the neutrino luminosities was reported 
also in Ref. ^8\. These features in neutrino luminosities are quantitatively the same for 
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Hyp-EOS and Shen-EOS models, and hence, it would be difficult to extract information 
of the NS matter only from the neutrino signal. Soon after the BH formation for 
Hyp-EOS models, /i/r neutrino luminosity steeply decreases because high temperature 
regions are swallowed into the BH, while luminosities of electron neutrinos and anti 
neutrinos decrease only gradually because these neutrinos are emitted via charged- 
current processes from the massive accretion torus. 

The anti neutrino luminosity for the long-lived HMNS for Shen-EOS models is 
Lp ~ 2 - 3 X 10^'^erg/s with small time variability. It is greater than that from 
the protoneutron stars found after supernovae [19]. Averaged neutrino energy is 
ep ~ 20 — 30 MeV. The sensitivity of water- Cherenkov neutrino detectors such as 
Super-Kaminokande and future Hyper-Kamiokande (HK) have a good sensitivity for 
such high-energy neutrinos in particular for electron anti neutrinos. The detection 
number for electron neutrinos is approximately estimated by aATLp/lAnD'^ep) where 
a is the total cross section of the detector against target neutrinos, AT is the lifetime 
of the HMNS, and D is the distance to the HMNS. For one-Mton detector such as HK, 
the expected detection number is > 10 for D < 5 Mpc with AT ~ 2 - 3 s, neutrinos 
from the HMNS may be detected and its formation may be conffimed. Note that 
GWs from the HMNS will be simultaneously detected for such a close event, reinforcing 
the confirmation of the HMNS formation. For Hyp-EOS models, the detection of anti 
neutrino is less likely because of the short lifetime of the HMNSs. However, as discussed 
below, the gravitational waves could constrain the hyperonic or nucleonic EOS. 

Figure [3] plots the plus and cross mode (/i+,x) of GWs with / = |m| = 2 as a 
function of tret — Emerge whcrc tret IS the retarded time tret = t — D — 2Mlog(Z)/M), 
extracted from the metric in the local wave zone, i.e., at ~ 550 km. The waveforms 
are composed of the so-called chirp waveform, which is emitted when the BNS is in 
an inspiral motion (for tret ^ Emerge), and the merger waveform (for tret ^ Emerge)- The 
GW amplitude is |/i+,x| ^ 2 x 10^^^ for a source at a distance D = 100 Mpc with 
the direction perpendicular to the orbital plane. GWs from the inspiral phase (for 
tret ^ tmerge) agree Well with each other for the models with Hyp-EOS and Shen-EOS 
with the same mass. On the other hand, quasi-periodic GWs from the HMNS (for 
tret ^ tmerge) show differences. Quasi-periodic GWs is suddenly shut down at the BH 
formation for Hyp-EOS models, tret — tmerge ~ 11 ms for H135 and tret — tmerge 
for H15. This is because the HMNS collapses to a BH due to the softening of the EOS 
before relaxing to a stationary spheroid. The quasi normal mode is excited after the 
BH formation, which is imprinted in the gravitational wave spectra at ~ 7 kHz (see 
Fig. H^b)). For Shen-EOS models, the long-lived HMNSs emit the gravitational waves 
whose amplitude gradually decreases indicating that they approach to the axi-symmetric 
quasi- stationary state. 

The characteristic GW frequency, few, increases with time for Hyp-EOS models. 
These facts are well imprinted in the effective amplitude (see Fig. IH^a)) defined by 
hes = 0.4/|/;,(/)| where h{f) is the Fourier transform of h+—ihx and the factor 0.4 comes 
from taking average in terms of direction to the source and rotational axis of the HMNS. 
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The peak amplitudes of hes{f) in Hyp-EOS models are smaller than those in Shen-EOS 
models due to a shorter lifetime of the HMNS. The prominent peak in the GW spectrum 
for Hyp-EOS models is broadened because of the shift of the characteristic frequency. 
This frequency shift is well imprinted in the evolution of the characteristic frequency 
/gw in Fig. H (b), where few = dcpcw /dt /2n with the Weyl scalar = Ae"*'^Gw_ Pqj~ 
Shen-EOS models, the central value of the frequency remains almost constant after the 
merger (tj-et ^ Emerge) and it corresponds to the peak frequency of gravitational wave 
spectrum in Fig.|l](a). By contrast, Hyp-EOS models exhibit a gradual increase in /gw- 
Specifically, for HI 35, /gw grows from ^ 2 kHz at tret ~ Emerge = 2 ms to ~ 2.5 kHz at 
tret — Emerge = 10 ms. For H15, it iucreases from ^ 2 kHz at tret — Emerge = 2 to ~ 3.4 
kHz at tret ~ ^mcrgc = 2.5 ms though the lifetime of HMNS is short. This feature in the 
frequency shift is explained as follows. For Shen-EOS models, the gravitational waves 
carry away the angular momentum and then the HMNS slightly contracts. For Hyp- 
EOS models, the hyperon fraction Xa increases with the contraction of the HMNS, 
resulting in the relative reduction of the pressure. This leads further contraction of 
the HMNS. Recent studies suggest that /gw is associated with the f-mode which is 
approximately proportional to ^/MnsT^ns with Mns and -Rns being the mass and 
radius of HMNS [201 EI]- Applying this result to our case indicates that Hyp-EOS 
models exhibit the frequency shift due to the run-away contraction, i.e., decreasing the 
radius of the HMNSs, and the Shen-EOS models emit the gravitational waves with 
the coherent frequency due to the slight contraction. In other words, the frequency 
/gw for Shen-EOS models evolves on a timescale of gravitational radiation reaction. 
The effective amphtude 1 - 4 x 10~^^ for D = 100 Mpc suggests that for a specially- 
designed version of advanced GW detectors such as broadband LIGO, which has a good 
sensitivity for a high-frequency band, GWs from the HMNS oscillations may be detected 
with S/N = 5 if < 20 Mpc or the source is located in an optimistic direction [6]. 

4. Summary 

We have reported the results of numerical-relativity simulation performed incorporating 
both a finite-temperature nucleonic EOS and hyperonic EOS as well as neutrino cooling 
effect. We showed that for both the nucleonic and hyperonic EOS, HMNS is the 
canonical outcome and BH is not promptly formed after the onset of the merger for 
the canonical value of the NS mass 1.35 Mq and for more massive case of 1.5M0. The 
primary reason is that the thermal pressure plays an important role for sustaining the 
HMNS. For the nucleonic EOS, the lifetime of the formed HMNS is much longer than 
its dynamical timescale, 3> 10ms, and will be determined by the timescale of neutrino 
cooling. For the hyperonic EOS, the HMNS subsequently collapse to a BH due to a 
softening of EOS as a consequence of increase of A hyperon. Neutrino luminosity of the 
HMNS was shown to be high as ~ 3 - 6 x lO^^erg/s. The effective amplitude of GWs 
is 1 - 4 X 10^^^ at /peak ~ 2 - 3.2 kHz for a source distance of 100 Mpc. 

We further found that the characteristic frequency of gravitational waves, /gw, 
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Figure 1. Maximum rest-mass density, maximum matter temperature, and maximum 
mass fraction of hyperons as functions of time for H135 (solid-red), S135 (dashed- 
green), H15 (short-dashed blue), and S15 (dotted- magenta). The vertical thin solid 
(short-dashed) line shows the time at which a BH is formed for H135 (H15). 



from the HMNS increases with time for the hyperonic EOS in contrast to the nucleonic 
EOS in which /gw approximately remains constant. This time-dependent characteristic 
frequency leads to the broaden gravitational wave spectrum for the hyperonic EOS. For 
the nucleonic EOS, the spectrum has a sharp peak around the less time-dependent 
frequency. This result suggests that the emergence of hyperons may be captured from 
the evolution of characteristic frequency and the peak width of the gravitational waves 
spectra. If the BNS merger happens at a relatively short distance or is located in an 
optimistic direction, such GWs may be detected and HMNS formation will be confirmed. 
In particular, to increase the chance of the detection of such GWs, we need for the 
detectors with high sensitivities at high frequencies of several kilohertz such as the 
Einstein Telescope [22]. 
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